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Abstract 
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the spin relaxation times for the Si-doped GaAs materials, obtained here 
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transitions. We also discuss the Mott transition within the context of the 
effective theory. 
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1 Introduction 



Both electronic and spin dynamics are major research areas in Condensed Matter 
Physics, specially for materials which induce strong correlating among electrons. 
In such systems, the interplay between several effective interactions and physical 
degrees of freedom is believed to lead to interesting effects, such as the colossal 
magnetoresistance and the high Tc superconductivity. The microscopic structure of 
these materials is complex: generically, these are compounds grew on some substrate 
with random impurities, yielding a crystalline structure with some amount of disor- 
der. Moreover, the low density of the free carriers renders often their kinetic energy 
comparable to the Coulomb energy, pushing the system close to the Hubbard-Mott 
transition {i.e., metal-insulator transition driven by correlations). In this regime, 
the approximations leading to the usual band theory results become invalid, and it 
is assumed that the electronic system is not a Fermi liquid. 

The Hubbard-Mott transition has been the subject of much attention in Con- 
densed Matter Physics, often studied in reference to the Hubbard Model. Actually, 
some authors recognize two types of Mott transitions [Tj: the so-called f/-Mott 
transition, in which at a given filling, a variation of the interaction leads to a metal- 
insulator transition, and the (5-Mott transition, for which the interaction is constant 
and the doping variable, leading to a commensurate to a incommensurate transition. 
The physics of the transition still remains not well-understood in d = 2 and d = 3 
spatial dimensions. Moreover, most of the studies used mean-field techniques, in 
particular, the Dynamical Mean Field Theory [2]. Some more detailed information 
of the transition, such as its universality class or critical exponents, still remain to 
be fully established. Recently, it has been claimed that the Mott transition is of the 
gas-hquid type, and belongs to the Ising Universahty class. [3] [1]. One may argue, 
however, that at least for some systems of interest, the separation between f/-Mott 
and 5-Mott transitions is somewhat artificial, since, at low enough density values, 
changing the number of free particles {i.e., doping) changes simultaneously the in- 
teraction and the commensurate-incommensurate character of the 6 Mott transition. 

On the other hand, the Mott transition has recently attracted interest in a rather 
different arena, namely, the non-equilibrium magnetic properties of spin polarized 
electrons. In a remarkable experiment, Kikkawa and Awschalom|5] have found that 
the spin lifetime in Si-doped GaAs semiconductors can be enlarged by up to three 
orders of magnitude by the introduction of ra-doping of the order of n = lO^'^cm'^. 
Among the early mechanisms that have been used to study the spin decay in semi- 
conductors are the ones proposed by Elliot-Yafet (E-Y)|6] D'yakonov-Perel (D- 
P)^ , Bir-Aronov-Pikus |8j, the relaxation due to nuclear spins, the Dyaloshinskii- 
Moriya (D-M) interaction [TU] [TT] and others. Recently, some calculations have been 
repeated for the E-Y mechanism, in the low temperature regime [T^, and for the 
D-P mechanism [13] but they disagree with the experimental results by nearly three 
orders of magnitude at doping near to n = 10^^ cm~^. Moreover, it has been claimed 
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that the enhancement of spin hfetime can be understood within the D-M picture |14j . 
This claim has been criticized by Gorkov, who has pointed out that the enhanced 
hfetimes take place at doping levels near to the metal-insulator (Mott) transition 
|15j . In fact, while all the preceding mechanisms and calculations are perturbative, 
the Mott transition is a strongly correlated effect, thus at first glance, the Hiller- 
London calculation in ref [H] seems misleading as was pointed out in ref [15]. At 
present, however, we are not aware of any calculations of the Mott transition for 
systems like GaAs. In particular, it is unknown how the crystalline symmetries and 
the band structure could affect the Mott transition in semiconductor systems. The 
answer to these questions cannot be given by performing band structure calculations, 
due to the limitations imposed by the strong correlations. 

To summarize: on the one hand, there is not as yet a consensus about the nature of 
the spin relaxation mechanism in n-doped semiconductors: several different models 
are all based in band structure and dilute gas scattering approximations, which 
may not be valid for a phenomenon involving strong correlations. On the other 
hand, several characteristics of the Mott transition are still not well understood, 
and its eventual relevance to the spin relaxation dynamics it is still an open issue . 
A first step towards the understanding of both problems would be the formulation 
of a unified framework for describing systems of electrons including both the band 
structure information and the strong correlations effects. The main purpose of this 
paper is to propose such a framework, given by an effective field theory that can 
be applied for crystalline condensed matter systems of the type discussed above. 
Effective field theories (EFTs) are a general, systematic and powerful framework 
to describe the low energy (long distance) dynamics of systems with one or more 
energy scales. These include and unify some earlier approaches, such as the Landau- 
Ginzburg theories and the Fermi liquid theories, among others [16]. They provide 
a systematic tool to analyze the different interactions and their relative importance 
through the use of the Renormalization Group (RG) ideas and techniques. Thus, 
this framework is well suited for the crystalline systems in the regime discussed 
above, where a naive calculation shows an energy scale of about E ~ h/lOOps ~ 
6.5 10~^ mev at doping levels of n ~ 10^^ cm~^, while the typical Fermi energy is of 
the order of Ep ~ 2 mev at this doping, and the gap energy is of the order of some ev. 
Eg 3 ev. As a test bench for the EFT, we discuss the Mott transition for samples 
of S'i-doped GaAs, and show that it can be used as a useful tool for discussing some 
properties of the spin relaxation dynamics processes, although we shall not address 
those processes directly, given that these are non-equilibrium process and therefore, 
outside the applicability range of any field theory. 

This paper is organized as follows: in the section 2 we write down the general 
EFT for electrons in crystalline systems. In section 3 we consider the inclusion of 
low disorder in the EFTs using standard techniques. Section 4 is devoted to the 
formulation and discussion of the EFT for the case of S'z-doped GaAs materials, 
including some perturbative Renormalization Group calculations. In Section 5, we 
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discuss the strong coupling regime of the EFT and its relation to the Mott transi- 
tion, the scenario for the appearance of a d-wave density order parameter and their 
relevance to the systems considered here. Finally, we give some Conclusions. 

2 Effective Field Theories for Lattice Systems 

Condensed matter systems are quite often hard to solve, due to their large number of 
degrees of freedom and to the strong interactions that arise among them. However, 
the description of many systems at a given energy scale becomes simpler in terms of 
a few representative degrees of freedom. This observation, which is a consequence 
of the Renormalization Group, is the key idea behind the formulation of a given 
system in terms of Effective Field Theories (EFTs) [15] [H] . The general scheme 
for applying the method of EFTs to a given system starts by conveniently choosing 
the characteristic degrees of freedom that describe the system at a given scale. In 
general, this choice is suggested by the phenomenology or the experimental data. 
Then, one proceeds to write down the most general local action which is consistent 
with the symmetries of the system, in terms of (second quantized) fields (or 'order 
parameters') that represent those degrees of freedom. However, one only retains 
a few of all the possible terms in the action, by analyzing the scaling behavior of 
each, i.e., only those terms that do not decrease in the limit of low energies or long 
distances (relevant and marginal terms). The so-obtained action is then called the 
(Wilsonian) Effective Active, and allows for a complete description of the system 
in the limit of low excitation energies, i.e., it yields its ground state and low-lying 
excitations. 

The systems we will consider here are non-relativistic electrons moving in crys- 
talline backgrounds. The basic field is taken to be non-relativistic fermionic (spin 
1/2) field. The basic symmetries of the system of electrons moving in a given back- 
ground are: i) the electron number, ii) the discrete lattice symmetries, iii) the spin 
SU{2) and iv) the U{1) gauge symmetry of electric charge. We would like to antic- 
ipate here that we will consider relativistic corrections to the standard Fermi liquid 
theory by allowing for spin-orbit couplings. Anticipating some of the discussion, we 
would like to make the following remark: in some crystalline systems, the full SU (2) 
symmetry group could he reduced to subgroups of it compatible with the symmetries 
of the point group of the crystal. 

The low-energy properties of the above electronic systems are described in terms 
of the dynamics of the Fermi surface, which is written in terms of the basic fermionic 
fields in the EFT formulation. The ground state of the non-interacting theory consist 
of a filled up Fermi sea with a Fermi surface. The elementary low-lying excitations of 
the system are particle-hole fluctuations above the Fermi surface. The geometrical 
shape and properties of the Fermi surface are important characteristics which serve 
for characterizing the system in a universal manner. For spherical (or convex) Fermi 
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surfaces, the EFT in three spatial dimensions yields the familiar Landau theory of 
Fermi liquids [TB],[Trj. However, there are other cases in which the interaction of 
the electrons and the crystalline structure can modify this simple picture, and it 
becomes necessary to take into account the electronic band structure. 

In order to illustrate some of the basic ideas of the method for constructing EFTs, 
let us first consider the following action: 

So = J dtd?x i)\x,t)[idt - H]i){x,t) , (2.1) 

where ilj{x,t) is a fermionic field and H the Hamiltonian of the system, which plays 
an important role in our non-relativistic systems. The field can be expanded in a 
Fock space basis as follows: 

^P{x,t) = J2 Mx,t) Cn , (2.2) 

n 

where c„ (cj^) are the fermionic annihilation ( creation) Fock space operators, re- 
spectively, and (j)n{x,t) are the corresponding wave functions. The Euler-Lagrange 
equations of motion are then Schrodinger equations: 

[idt - H](Pn{x,t) = (2.3) 



For time-independent Hamiltonians and translation-invariant crystalline back- 
grounds, the solutions to (12. 3p may be chosen in the Bloch form: 

0„(x) = e'^'-^-Hr) , (2.4) 

where k is the wave-vector, which belongs to the first Brillouin zone. The spatial 
periodicity then gives rise to the band structure, defined by: 

HM^) = e„0„(k) . (2.5) 

A functional integral of the fields that describes the system is given by the partition 
function: 

D'ip^Dip exp (^i^ j dt J d^k iplik, t) [idt - e„(k) - /x] ?/^„(k, t) 

' (2.6) 

where /i is the chemical potential, and the Fi denotes the first Brillouin zone. 

We would now like to construct the EFT corresponding to (12. 6p . We consider 
a characteristic energy scale -Ea, and a corresponding set of wave momenta at the 
Fermi surface k such that |k| < A. Here A is a chosen cut-off scale that serves to 
divide the fields into their high and low frequency components: ip = ipi + '4^h- More 
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precisely, 'i/'L(k, t) depends only of those momenta k satisfying |k| < |A|. The low- 
energy EFT is obtained by formally integrating out the high frequency components 
in dri]): 

Di){D^L e*^^(^^^^) , (2.7) 



where 

^ j Di,%Di,H exp(^5o(V'LV'L,V'LV'/f)) • (2.8) 

is known as the low energy or Wilsonian effective action. It can be expanded in 
the space of all theories in the (infinite) basis of all local operators On allowed by 
the symmetries: 

5a = I dt I d^k On . (2.9) 

•J -J n 

Notice that this effective action may contain free and (local) interacting terms among 
the chosen effective fields. For the case at hand, given by (12.71) we obtain: 

Z = j Di,{Di,Lexvi^Y. j '^^ d'ki,ln{^,t)[idt - e„(kL)-^]V^,„(k,t) 

(2^10) 

where Fa is not the entire Brillouin zone, but only the region near its center satisfying 
|k| < |A|. Therefore, all momenta appearing from now on are small, and it is possible 
to expand any function of momenta, like the energy, in a power series as follows: 

^^n(k) = Oii ki -\- Oiij kikj + Pijk kikjkk + fiijk kikjkj^cji + . . . , (2.11) 

where the Greek quantities denote coefficients and the Latin indices ran among the 
three spatial directions {i = 1, 2, 3). The different coefficients appearing in (12. lip are 
actually spatial tensor structures, which can be classified by the spatial symmetric^ 
of the system. In fact, it is well known that near the center of the Brillouin zone, the 
electronic wave functions transform as the irreducible representations of the point 
group of the crystal. 

In the following, we will focus on the EFT only, and therefore will suppress the 
subindices A and L in the effective action and effective fields, respectively. Therefore, 
for low energy and large distances the free part of the effective action may, at the 
n — th band, be written as: 



5^") = JdtJ d^ki;l{k,t) 



MKt) , (2.12) 



^Note that only third order tensors can couple the spin operator, because these may transform 
as a pseudo-vector in one index and as a vector in the remaining two, giving rise to a scalar when 
contracted with the spin pseudo-vector 
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where runs over all the different representations of the point group and the en- 
ergy operators E have to be expressed in terms of the basis functions within each 
irreducible representation. This generalizes the case considered above for more com- 
plex Fermi surfaces, that may display less familiar dispersion relations and several 
components. The elementary excitations are particle-hole pairs above the Fermi sur- 
face, which is defined by the equation e„(k) = ej?, whereas ei? is the Fermi energy. 
For spherical Fermi surfaces, these elementary excitations are Landau quasiparticles 
(dressed electrons) as has been pointed out in ^17] . The fact that the Landau theory 
of Fermi liquids is an EFT has been discussed in detail in [16] . In order to separate 
the different scaling components of the momenta, we consider the decomposition 
k = ko + p, whereas ko is a vector lying on the Fermi surface and p is a vector 
orthogonal to it. The scaling down of the energy E — >■ sE, with s — > a small di- 
mensionless parameter, induces a corresponding scaling on the different components 
of the momenta such that |k|o |k|o and |p| -sIpI- The energy operators are 
constructed accordingly: 

Er^{ki,kj,kk,. . .) = Er^{koi,koj,kok, . . .) + — M^^oj^^ofev^ ^^2.13) 
There is a set of different Fermi velocities given by the last term of the R.H.S of 

dan. 

We now turn to the interactions. Note first that the aim of the EFT approach 
is not to deduce the form of the interaction terms in the effective theory ( i.e., 
the interactions among the effective degrees of freedom) by direct transformation 
of an underlying microscopic interaction, but rather to make a judicious ansatz of 
the possible interactions in terms of the effective field symmetries (for a detailed 
discussion see [16]). At zero temperature, fi = ep and the dynamics of the low- lying 
excitations {i.e., the EFT fields defined as discussed above) above the ground state 
is given by: 



So = J dtj d^k^Pl{k,t) 



^ dEr^{koi,koj,kok, ■ ■ ■> 
lOt - 2^A« — Pi 



V'n(k,t) 



dEr^jko i, kpj , ko 

dki ^ 

(2.14) 

where the values of the dimensionless parameters Aq, are not predicted by the EFT, 
but are expected to be of order 1 ('naturality' [16j). If one tries to go beyond the 
scope of the EFT formulation, and insist in matching it with an underlying short- 
distances theory, then one may conclude that they depend on the details of the band 
structure and may be exactly zero only if there is an additional symmetry, or very 
small if the band is very far from the energy scale one is considering. 

The interaction terms may also be included in the EFT in a similar fashion. For 
convenience, we consider here the fields as functions of frequency, rather than the 
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time coordinate (using the obvious notation (1) = (cjijki)): 



(2.15) 



with the integration measure d^i = duJid^kiduj2(Pk2duj^d^k2duJid'^ki and the short 
hand notation: 5k5^ = 27r5(/ci + ^2 — /cs — k4)5{uji + UJ2 — oj^ — u^) where 5 is a 
deha-function enforcing momentum conservation mod 271 /a (a is the cubic lattice 
spacing), which allows for Umklapp process as well: 



ki + ka 



Q 







(2.16) 



where Q = m(7r/a, 7r/a, 7r/a) and m is an integer . Including higher interaction 
terms is not necessary due to the more and more irrelevant character of them jl6] . 

The complete EFT action is thus obtained from f l2.14p and (12.151) . and it is 
formally given by the following partition function: 



Dip^Dipexp \ i / dtd^kiljl{k,t) 



iOt - K eh: Pi 



dki 



V'n,(k,t) + 



/ c?/i[V^t^(l)^t^(2)V^„3(3)V^„,(4)\/(ki,,k2,k3,k4)]5^5. 



(2.17) 

At this point, we emphasize that the field operators ?/'„^(k, t) do not represent the 
original (free) electrons, but rather the effective (interacting) ones. 



3 Effective Field Theory for Impurity Systems 

In this section we would like to discuss the effects of impurities and doping in the 
systems we have considered in the previous section. The approach we take here is a 
standard one, based on the quantum mechanical established procedures which are 
promoted onto the field themselves by considering them as a linear superposition 
of the second-quantized operators acting on a Fock space, with amplitudes given 
by the modified wave functions. Quite often, impurity effects lead to interesting 
phenomena, such as in the case of high Tc superconductors in underdoped cuprates. 

The inclusion of impurities makes necessary to reformulate the EFT, since they 
break translation invariance, so that Bloch wave functions are not solution of the 
Schrodinger equation of the crystal. However, since in the EFT formulation we may 
to choose the relevant degrees of freedom, we consider the 'impurity orbitals' as the 
new basis for the wave functions appearing in the field definitions. They have the 
important property that can be extended to the mesoscopic domain. In fact, as is 
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well known, for a single impurity the Wannier theory of shallow donors shows that 
the hydrogen -like orbitals of typical impurities in semiconductors have an effective 
Bohr radius of about 50 Aup to 100 A. In order to show how this mesoscopic scale 
arises and how it matches the EFT framework, we would like to recall briefly the 
Wannier theory of shallow impurities [21] [28] [20]. Let us consider an impurity 
characterized by the potential SV, within a perfect crystal (quantities refer to it 
with a subfix). A solution of the time dependent Schrodinger equation 

[Ho + 6V{r)] 0(r, t) = thdt(j){r, t) (3.1) 

may be expanded in the Wannier representation 

0(r,t) = ^/„(ri,t)a„(r-r,) (3.2) 

n,l 

where a„(r — ri) are Wannier functions {i.e., discrete Fourier transforms of the Bloch 
functions) which are solutions of the pure crystal electronic problem, ri is the real 
space position of the atom, and fn{fht) play the role of enveloping functions. After 
a standard transformation, this leads to the equation (for details see [28]): 

[eni-iVi) + 6V{r)ni,n'i'] fnirht) = ihdtfnirht), (3.3) 

where e„(/c) are the band energies of the perfect crystal, V; is the discrete gra- 
dient with respect to the impurity site coordinates, and SV{r)ni^n'i' is the matrix 
element of the impurity potential between different Wannier functions. Now 6V{r) 
may be expanded around r; , under the assumption of a slowly variating potential, 
ao|V5l^(r)| -C 6V{r) and this equation becomes 

[enHVi)+SViri)]Ur,,t) = th dtUr,,t) (3.4) 

Note that equation (13. 4p does not involve the Hamiltonian of the perfect crystal 
Hq. Knowing e(k) we can solve for /„(r,t) in a continuous approximation. Since 
e(k) contains the symmetries of the lattice, this will also occur with /„, provided 
SV is small. If, as it occurs in a semiconductor, the band energy is still parabolic 
and the potential corresponding to a single 'shallow impurity' then the equation 
becomes the Schrodinger equation of an Hydrogen system, except that the Coulomb 
potential is replaced by an screened effective Coulomb potential, whose dielectric 
constant is modified by the presence of the host crystal Vg = e/Aner. The change in 
the dielectric constant gives raise to new spacial scale for the wave functions, and 
the effective Bohr radius becomes of the order of 100 A, well within the mesoscopic 
scale. 

To write down a field theory, let us consider an 'arrange of single impurities'. 
The degrees of freedom we chose are associated to the impurity level, so that the 
enveloping wave functions will be a suitable basis, z.e., the field operators are: 

^(r,t) = 5^c„/„(r,t) (3.5) 
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where c„ and cj^ are the annihilation and creation operators of electrons in the effec- 
tive bands whose wave functions are /„. An equation like (13 .41) may be considered 
as the Euler-Lagrange equation corresponding to the following action: 

So = J dtd?k V-Kk, t) [tdt - e„(k,) - 6V{-k) - /i] Vn(k, t) (3.6) 

where (5V^(x) = (5V^(x + R) plays the role of a 'lattice potential', now defined over 
the mesoscopic domain, which is invariant under the discrete translation symmetry. 
Since the functions e„(k) are the exact energies of the pure crystal, they contain 
the microscopic symmetries that near k ~ are those of the point group. In this 
way a single impurity extend these symmetries over mesoscopic domain. In this way 
we have recovered the discrete symmetry, so the electrons can exchange discrete 
momentum with the 'impurity lattice'. Note that 6V transforms as the temporal 
component of a gauge field Aq: 

So = y^ [ dtd^k^lik, t) [iD, - e„,(ki) - /i] V^„(k, t) (3.7) 

where Dq = dt — iAq . In real samples, however, impurities are not localized in a 
perfect spatial arrange, but are randomly distributed in a complex way that depends 
on the impurities, on the host crystal and on the growing process, so that the 
behavior of the EFT could depend on some of this data. Here we shall consider only 
the case of substitutional shallow impurities, such as those that appear in samples 
grown by molecular beam epitaxy (MBE). For this case, we have a 'low degree of 
disorder'; since the impurities are placed in some of the sites of the crystal lattice, 
the Hamiltonian of the problem can be modeled by: 

H, = Ho + g{Q) 5y(xi) (3.8) 

with g{(i) = 1 if Xi = and 5'(Ci) = if Xi 7^ ^i, and Q denote the impurity position. 
Now we have to solve the eigenvalue problem: 

HM=eM (3.9) 

where g acts as a parameter that takes values in a stochastic fashion with some 
definite probability distribution P{g). It follows that the set Hg may be considered 
as a ensemble of Hamiltonians in the random matrix theory (RMT) sense, with 
a probability distribution fi{H) induced by P{g)- As is well know in RMT, the 
statistical properties of the ensemble are fixed by the symmetry properties of the 
Hamiltonian ensemble. Moreover, the classical Wigner ensembles are defined by 
invariance requirements |TH]: given a particular member H in the ensemble, all 
the others are generated by similarity transformations on the Hilbert spaces. This 
postulate ensures that there is no preferred basis on the Hilbert space. Following 
this line of thought, we assume that the symmetry of the random matrix ensemble 
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is the same as that of the EFT. Physically, this means that we shall assume that the 
impurities form an array, such that the symmetries of the crystalline structure are 
extended onto the mesoscopic domain. A similar mechanism has been argued in [12] 
for the case of cuprate superconductors, in the sense that the d-wave order parameter 
survives to the presence of disorder. We shall discuss briefly the applicabihty of this 
hypothesis to the case of S'i-doped GaAs in the next section. 

A more precise formulation of the previous ideas starts by considering equation 
for fixed g, which defines the action: 



Sog = j dtd^'kijiik, t) [idt - gAo - e„^(k) - /i] ^„(k, t) . (3.10) 

Following the same steps that yielded (12.111) . we obtain: 

en,g{k) = afki + afjkikj + Pfj^^ikjkk + ps^jj^kikjkkai (3.11) 

Since we have assumed that within each array the lattice symmetries survived the 
change of scale, the energy operators are still defined by the irreducible represen- 
tations of the point group of the crystal, so that the effect of the disorder appears 
throughout the coupling constants only. The partition function Z\ for the theory 
parametrized by a particular value of A is given by: 

To complete the definition of the EFT, we must average over the disorder parameter 
A. This can be done in at least two cases: i) quenched disorder, which for static or 
thermodynamics theories is often handled thorough the replica trick. Such is the 
case,e.g'., of the Sherrington-Kirpatrick model of spin glass. or the field theoretical 
calculations of the average density of states in the Anderson localization problem. [18] 
|22j . ii) dynamical theories of spin glasses avoid the use of the replica trick. In 
these theories, the generating functional is normalized [23] [21] and the quenched 
average can be done directly using the generating functional. We will present here 
an approach to averaging over the dynamical theory that is inspired in the 'static 
formalism' and follows closely the approach presented in [25] . 

Let us consider the following generating functional, for a fixed set {Aq, Ar„}: 

DipDip^exp{t J dt£k{Leff + ^\k,t)J{t) + ij{k,t)J{ty}), (3.13) 
with 

L,jf = ^\k,t) ^tdt + \oAo-J2>^aErAki,kj,..)'^^{k,t) (3.14) 
For fixed A, the Green's functions are defined by: 

Gx{t - t', k) = {T{ij\k, t)^(k, t')) , (3.15) 
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and the average over the disorder (denoted by double brackets) can be done using: 

(3.16) 



6J{t')6J{t) 

where G\{t — t') = J d^kG\{t — t',k) . We now apply the rephca trick, i.e., the 
identity 

InZ = linin^o 



n 



(3.17) 



The average Green's function ( denoted by G]^{t) with r = t — t') becomes 

1 /■ " 

G'k(r) = d.jtdMdim^_o- / \{Di;,Dil,ld\P{\)d^^''^ (3.18) 

^ a=l 

S[\] = J dtd'kil^l [^^t - AoAo - KEtJ\ + J^i + JHa (3.19) 

where we have used djtdjf = 5"^ / 5J{t')5J{t). We consider a standard Gaussian 
probability distribution for the coupling constants. 



P{\) 



1 



,-(A-Ao)2/(2^') 



2Txa 



(3.20) 



so that the integration over disorder is Gaussian, and can be done in a closed way: 

1 



G = djtdjt'limn^o- 



n 



a=l 
2 



II DiPaDiljl exp / dtd^'kijl \idt - AoAo - K^Et} + -Hi + JHa 



af3 



^^PtE / rfti«VA;2^1(l)V^a(l)i5;{^Jl)A^(2)V^,(2)E{^^(2)) 



exp 5^ / c?/i<(lX(2)V^„3(3)7/.„,(4)r(ki,k2,k3,k4)45. 



(3.21) 



where we have included the interaction term, since we assume that the Coulomb 
potential is independent of the disorder. In the second line of 03.211) . we have used 
the short-hand notation E-pg = Aq. Note that the impurity potential can mix with 
other operators. 

Thus, integrating out the random disorder we have introduced an effective inter- 
action between operators with different symmetries. Effective interactions arising 
from disorder were already present in early models, such Sherrington-Kirkpatrick 
model, and in the study of the classical dynamic of spin glasses [23] • In the problem 
at hand, they acquire a new meaning. In fact, it has been argued that in disordered 
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or complex condensate matter systems, different symmetries and order parameters 
may compete [2S]- In the problem we are considering, this behavior results directly 
from EFT: we have began with a set of definite degrees of freedom and given sym- 
metries but the process of averaging over the disorder yields a theory mixes those 
degrees of freedom and new induced interactions, which compete with the Coulomb 
interaction, arise among them. For the rest of the paper, we will consider the case of 
low disorder only,ie., a small so that the effective induced interactions are negligible 
in comparison to the Coulomb term. 

4 Effective Field Theories for systems of 52-doped 
GaAs compound materials 

In this section, we apply the framework of the preceding one to analyze the behav- 
ior of the S'i-doped GaAs samples employed in the spin relaxation experiments of 
reference [5]. We obtain the effective Hamiltonian and analyze the RG lowest order 
perturbative corrections to the coupling constants values. 

Let us first consider an array of Si impurities in GaAs with the symmetries of the 
substrate. This hypothesis is reasonable, because GaAs has zincblende crystalline 
structure and the silicon grown in fits in a diamond lattice (with the same structure 
of the zincblende, but with only one kind of atom in the basis). The samples of 
S'i-doped GaAs considered are grown by Molecular Beam Epitaxy, which implies 
that the impurities are almost all substitutional. Moreover, the single-impurity 
wave functions have sizes of about 100 A, and for doping near lO^^cm"^, these are 
separated by distances of about 200 A. In addition, some self-organized structures 
and three dimensional nanowires have been observed inside samples of Sz-doped 
GaAs [22] . Note that this hypothesis may or may not be valid for systems different 
to the one considered here. A general rule as to how to proceed cannot be formulated 
at this point and the specific properties of each system should be taken into account 
(as is often the case in the framework of EFTs). 

From the symmetry point of view, the zincblende belongs to the spatial group 
or ^43^ [28], which has 24 elements (it is the same as the point group of the 
tetrahedron, ). There are 8 C3 rotations about the diagonals (±1, ±1, ±1), 3 C2 
rotations about the axis (1, 0, 0), (0, 1, 0) and (0, 0, 1), and 6 C4 improper rotations 
about (1, 0, 0), (0, 1, 0) and (0,0,1). The group Tj, has also six reflections a with 
respect to the planes [1, 1, 0], [1, 1, 0], [1, 0,|l]], [Oil] and [Oil]. On the other hand, 
Silicon belongs to the spatial group 0\ (the diamond group), which is non-symorphic 
and has three glide planes defined by x = a/8, y = a/8, z = a/8 and screw axis 
along the directions x, y, z. 

li is well known that electrons in Si are located in impurity levels near the conduc- 
tion band (donor levels) (and the holes are near the top of the valence band, accep- 
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tor levels) There are six Fermi points or 'pockets' along the directions (±1, 0, 0)7r/a* 
,(0, ±1, 0)7r/a*,(0, 0, ±l)7r/a*, with energies: 



eik) = e, + ^(^ + ^ + ^) (4.1) 
2 \mi rrit rut J 

where mi and rrit are the effective masses in the longitudinal and transverse direc- 
tions. As discussed in section 2, in order to write down the EFT one considers the 
low-lying excitations above the ground state, which can be done (at lowest order in 
momentum) by linearizing the dispersion relation about the pocket points yielding 
three (six chiral) quasi one-dimensional gapless systems (Luttinger liquids, see Fig- 
ure 1). However, as we shall soon see, there are more allowed terms in the EFT 
besides these ones. 

In order to obtain all the terms, we apply the approach of section 2 and construct 
the energy operators Er^ ? which transform as the irreducible representations of the 
symmetry group of the system. In our case, we consider energy operators invariant 
under the symmetries of both Silicon and Ga — As,i.e., diamond and zincblende 
structure. But since is a subgroup of the diamond group, all the symmetries 
needed to construct the EFT are in fact contained in the group T^, which has the 
following characters and basis functions table: 





E 


3C2 




6a 


8Cs 


Basis Functions 


A, 


1 


1 


1 


1 


1 


kx ' ky . /c^ 




1 


1 


-1 


-1 


1 


kt{k1 — k1) + ky{k1 — k1) + k^{k1 — ky) 


E 


2 


2 








-1 


{{kx — ky), 2k^ — {k^ + ky)} 




3 


-1 


1 


-1 





{kx{ky — k^), ky{k^ — k^), kz{k^ — ky)} 


T2 


3 


-1 


-1 


1 





{kx, ky, kz} 



Note that only the class Ti transform in a pseudo-vector representation, so it must 
couple to the spin variable to produce a scalar term in the effective action. 



Using only the equations (12.130 . the character table and neglecting four order 
terms {i.e., those arising from the A2 class), we obtain the effective Hamiltonian 
Hq^^ = Hp + Hgo, where the second term is the spin-orbit Hamiltonian: 



Hr, 



H. 



XqAq + Xi\Px + Py + Pz 



[X2{gxa-P) + A3(gs/Q.p) + X^igza-P) 



kox PyPz ~l~ kgy PxPz ~l~ kgz P xP y 



koxipl - pI)(Tx + koyipl - pD'^v + kozipl - pl)(y. 



(4.2) 
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where (Ajt, k = 1, . . . ,6) are coupling parameters, gx^± = (0, 1, ±1), gy^± = (1, 0, ±1), 
gz^± = (1, ±1, 0), which are vectors in the reciprocal space. Note that these vectors 
appear also in the pseudo-potential theory for GaAs ^28j . In (14.21) one may consider 
the terms linear in momenta as subsystems of quasi one- dimensional Luttinger liq- 
uids. There are also two parity-breaking terms: the spin-orbit part has the form of 
the Dyakonov-Perel Hamiltonian (also called Dressehauss term) [9j. Note that term 
of the form {k^y — kl^)pxCrx, and the ones obtained from it by cyclic permutations 
could also be considered in Hgo- However, they do not contribute to Hgo for the 
following reasoning: a parity transformation acting on the low-lying states of the 
theory not only transforms the components of the electron momentum p, but also 
exchanges two of the corresponding pocket ground states of Figure 1. The parity 
transformation reflecting the axis x onto minus itself can be seen as an exchange of 
the axis y with z, while keeping x fixed. Therefore, the term [k'^y — k^^) changes 
sign, as well as ax, while px does not. However, this term is identically zero since 
k^y = k1^. Moreover, one could also consider a chiral component of each Luttinger 
liquid {i.e., one of the six pockets of Figure 1), and after adoption of a convenient 
redefinition of the parity transformation on this system conclude that these terms 
neither arise in this case. 

The field theory corresponding to the Hamiltonian (14. 2 p is defined by its partition 
function, given by: 

Z= j D^[DtljLe'^ . (4.3) 



The free part of the action is given by 

^0 = E r {i^K^^Pc.) h - ^0^0 - XPa + Oip')] U^,p^)} (4.4) 

where the index i runs over the six Fermi points in the different spatial directions 
{i.e., i = Rx, Lx, Ry, Ly, L^, Rz, where R, L mean left and right Fermi points, respec- 
tively), and a is an index that indicates the component of the momentum along the 
six different axis defined by the directions x, y, z, Qx, Qy, Qz- By abuse of notation, we 
take a = x,y,z,gx,gy,gz- Note that the two indices i and a) are not independent; 
for example, ii i = Rx then a must take on of the values given by x,gy or gz- In 
Figure ([T]) we indicate each of the six possible Luttinger ground states ('pockets') 
and a few possible intervalley processes produced by the effective Hamiltonian (14.21) . 



The general form of the interaction term is 
= 2!2! E ' ^^^l(^)^l*^2)^^3(3)V'i4(4)^ii,i2,i3,i4(^oi,/co2,/co3,/^04)45^ 
where 

du 



"'^^ dpa,i f duji r^=^ dp% f du2 dp^^^^ f dco. 



Pi=-A 27r J 2tt Jp^=^A 27r J 2tt Jp^=-A 2tt J 2-n 
Sk = S[uii{koi +pi) + i^i2{ko2 +P2) - i^i3{ko3 +P3) - z/i4(A;o4 +P4)) , (4.5) 
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Figure 1: Ground state structure of the EFT in terms of Luttinger liquid 'pockets'. 

with Ui = ± when i = R,L and = 6{u!i + UJ2+UJ3 + UJ4). It is not difficuh to show 
( for example, in a real space nearest-neighbor model, for details see [I7j) that the 
interaction momenta-dependent couplings are of the form 

V+{kik2k3, ki) = Uo sm( — - — ) sm( — - — ) cos I 1 (4.6) 

with k = ko+p, and where Uq is an ultra local (on-site) repulsive Coulomb potential. 
Two main possibilities exist for the scaling behavior of the coupling functions in 
(14.51) . If ki and kj arise from opposite points in the Fermi surface, 

V ~ sin( ^' ~^' ) ~ sin(7r/2 + 0{p)) (4.7) 

On the other hand, if ki and kj arise from the same side of the Fermi surface 
(Umklapp process) 

Vumklapp = {Pl - P2){PZ " Pi) ~ 0{p^) (4.8) 

We are now ready to study the behavior of each term in the effective action for small 
energies and momenta (relative to the Fermi surface) using the RG techniques. The 
change in scale is given by: 

Uj' = su 
p' = sp 

x' = s-^x (4.9) 
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where s is a small parameter such that s — 0. The fields scale according to: 

^\u[p[) = s-'/'^icOi^pi) , (4.10) 

so that the terms linear in p (and the frequency terms) are marginal, while the terms 
0{p'^) are irrelevant. 

The scaling of the interaction terms (which are free of ko) can be deduced from 
the scaling properties of the delta-functions: 

Kp)^5{v'/s) = s5{p') (4.11) 

5{ij) 5{u'/s) = s6{uj')] (4.12) 

Then the inter-valley processes (those in which fc, and kj come from different Fermi 
points) are marginal. The Umklapp process are suppressed by a term 0(p^),z.e. they 
scale as and are, therefore, irrelevant. The scaling of the gauge potential depends 
on the explicit form of Aq. For shallow impurity donors, it scales as ~ l/\r — r'\ 
(Hidrogen-like Coulomb potential) so that it goes as in three dimensions, and as 
in two. Hence, they are also irrelevanlH. 

Hence, the complete EFT is given by the following action: 

i,a 

+ f rf/^^i(l)^l(2)V^^.3(3)^M(4)2vrVc. (4.13) 



Here ^ is a coupling constant (Umklapp processes are excluded). The simple pic- 
ture that emerges is displayed in figure ([T]): the low-energy dynamics of the three- 
dimensional Fermi surface can be described in terms of six one- dimensional Luttinger 
systems. 

Since the quartic terms in the EFT action (14.141) are marginal, it is necessary to 
consider their quantum corrections, i.e., to include one-loop contributions in the RG 
calculations. This is done by writing down all the Feynman diagrams up to second 
order in the interaction coupling constants and integrating out all momenta in a 
shell of width AA. For one-dimensional systems, this integration is made about the 
four Fermi lines (labeled a, b, c and d in Fig(l2l)) located at distance A from the Fermi 
points with width AA. In the regime of small momenta, the Coulomb potential is 
usually taken to be constant, and all the external frequencies and momenta are set 
to zero, obtaining 

AV = ^ [ diu [ dp- ^ 



4vr2 7 '^[tuj-E{p)][tuj-E{p)] 



^Ultralocal potentials of the form d6{x — x')/dx that also break parity are be marginal. 
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A 



Figure 2: Regions of momentum space around each Fermi point, defined by a, b, c, 
d, being integrated out. They he a distance A from the Fermi points L and R and 
have a with AA 



47]-2 



dp- 



1 



1 , 



AA [iuj - E{p)] [iuj + E 
1 



dp- 



(4.14) 



471^2 J '^[tuJ-E{p)][-tuj-E{p)] 
where we have used the symmetries E{p) = E{—p) and E{p±it) = E{p) (see figure 



Evaluation of the integrations leads to the Beta function j3 = AdV/ dA: 



PiV) 



2ti 



27r 







(4.15) 



This result implies the existence of a line of continuous fixed points ( see, e.g. 
|17j). The emerging scale invariant states are identified with those of a Luttinger 
liquid. It follows from the above calculation that the Luttinger liquid phase is 
stable against quantum corrections. This perturbative analysis makes sense if the 
low-energy effective dynamics picture of figure ([T]) is stable under perturbations, 
which is a hypothesis we have implicitly made. 



5 Strong coupling regime and the Mott transition 



As we have discussed in the previous section, the EFT of the three-dimensional 
crystalline structure of S'i-doped GaAs is given by a system of six coupled one- 
dimensional Luttinger liquids, which include Umkplapp processes. In order to gain 
some insight about the behavior of the system, we will neglect for the moment 
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the coupling among the Luttinger hquids, so that we can focus on reduced one- 
dimensional dynamics. Each Luttinger system is exactly solved [22] and from their 
exact solution, it is known that the excitations of the system include charge density 
waves. From the RG point of view, as one moves along the line of fixed points, 
the scaling dimensions of the various fields change continuously away from the free 
field fixed point values. Simultaneously, the Umklapp processes that were initially 
suppressed by a factor 0{p^) change their scaling behavior, firstly becoming marginal 
and later on, relevant operators. In this new scenario, their effects should be also 
considered in the low-energy description of the system. 

In generic conditions, the Umklapp process are short-distance, high-energy effects 
which are not considered in the EFT, since the scale of the involved momenta is of 
the order of vr/a (a is the microscopic lattice spacing). In fact, the low-energy 
physics depends on short distance dynamics only through relevant and marginal 
couplings, and possibly on a few irrelevant terms, if one measures small enough 
effects [IGj. However, for the class of systems we are considering in this paper, the 
low disorder impurities extend the microscopic lattice symmetries to a mesoscopic 
domain, implying that Umklapp processes are also promoted to mesoscopic scales. 
The transferred momentum between lattice and electrons at these scales is now of 
the order of vr/a*, where a* is identified in the physical system with the effective 
Bohr radius of the impurities (a* ~ lOOA). Note that this mechanism is not fine 
tuned since a* depends on the screening and therefore, on the density. There will 
be a range of values of a*, given by a corresponding range of the doping density 
n, for which Umklapp processes should be included in the low energy theory. The 
appearance of Umklapp process in this systems is also interesting because they have 
been put in correspondence with the Mott transition. Their effects can be naturally 
included in the EFT after bosonization of each Luttinger liquid. 

One dimensional systems may be described in terms of density and spin fluctua- 
tions (see, e.g., [30] and 



H — Hp + H„ + H scatter + H umklapp (5-1) 

with 



Hj, = — I dx 

271 



Humklapp = ^2nl\'2 J ^^CO^[^<^P + "^^i 



(5.2) 
(5.3) 



Hscatter-- ^^^^^ 



^ j dxcos[V8(l)^] . (5.4) 



Here u = p or a denote the charge and spin labels, and the bosonic fields (j)p and 
0O- describe the charge and spin fluctuations, respectively. The parameters Kp, K„, 
Up, Ua and are the only necessary data to describe the long distance physics. 
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The parameter 5 = Akp — 27r/a is measure of the distance away from half filhng, 
and n/2a is the Fermi vector at half filling. Therefore, the doping d is related to 6 
through d = a6/2n. The RG equations in these variables yield the system: 



^ = -\y!K'M6{l)a) (5.5) 

^ = (2 - 2K)y, (5.6) 

^ = -lyluKJMl)a) (5.7) 

f = m + ^yPMl)a) , (5.8) 

where = g^/{'Ku) and I describe the renormalization of the cutoff a{l) = aexpl, 



and Jj, are Bessel functions. For half filling, the RG equations coincide with those 
of Kosterlitz-Thouless, with a separatrix at K — 1 = g^/{2'Ku) between a regime 
for which is irrelevant and another one for which it becomes relevant. Given 
that vortex condensation or limiting surface creation are known to arise in systems 
for which the Kosterlitz-Thouless transition develops, one might expect that these 
physical processes could similarly appear in our case. 

The Mott (or more properly, the Hubbard- Mott) transition is a strongly correlated 
effect yielding a metal-insulator transition when the electron interactions become 
strong enough. Very frequently, it is studied within the realm of the Hubbard 
model. 

H = -Y^ t{clcj„ + c]^Ci„) + UY^ rii^rii^ , (5.9) 

where t is the hopping parameter and U the on-site Coulomb repulsion. The tran- 
sition itself is interpreted as arising from the competition between the kinetic term, 
which favors the delocalization of electrons and the interaction (repulsive) term, 
which favors localization. Since equations |5.2j . |5.3j . |5.4j . are the most general for 
one- dimensional fermionic systems, they describe the Hubbard model in one dimen- 
sion as well. If fact for weak coupling, it is possible to show that: 

UpKp = u^K^ = VF Up/Kp = VF + U/'K (5.10) 



uJK^ = vf-U/tx g3 = U. (5.11) 

Therefore, the interaction driving the Umkplapp processes is the on-site Coulomb 
repulsion. Let us consider the charge sector of the one dimensional theory near half 
filling. Parametrizing the deviation from this point through the chemical potential 
we can write the Hamiltonian as: 



H 



— I dx 

2-K 



UpKp{-KT{p 



Kp 



+ pl^Cos[V^<Pp]-lid, 
2Tia'^ 



(5.12) 
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After introduction of the new field, (j) = \f2(\)p it is possible to rewrite the Hamilto- 
nian (for details see [SU], [22]) as: 



H = J2'"^(^kCL,k - c[kCL,k) + -^Y1 ikCR,k + h.c 

k k 

^TmShm J2 2pMP-i-p) - fMp)P+{-p) + P+(p)p+(-p)) (5.13) 
p 

where Ps{p) denotes the Fourier transform of the fermionic density {s = ±), v = 
?/(cosh2^ + /isinh2^) and exp(2e) = l/i2Kp). For the value Kp = 1/2, the so- 
called Luther-Emery line, the two fermion species become non-interacting and the 
Hamiltonian can be diagonalized by a Bogoliubov transformation, given raise to two 
bands in the spectrum. 



El,r = ±Vivk)^ + ((?3/(27ra))2 (5.14) 

Therefore, the gap energy is given basically by the the Umklapp coupling constant. 
Away from this line, the interaction term gives rise to a renormalization of the gap. 
The energy values are: 

EL,R = ±y^{vky + {A/2y , (5.15) 

where A is the renormalized gap. For enough small doping, (near half filling), the 
chemical potential is located in the gap region and at zero (or low )temperature, 
the system becomes insulating. As soon as the chemical potential crosses the upper 
(or lower) band, a chiral current is established. This is a brief account of the Mott 
transition in one dimension. In two dimensions, further difficulties arise. However, 
several techniques have been developed to deal with strongly coupled fields for these 
systems. The analysis of the non-trivial fixed-points of the RG in this case is beyond 
the scope of the present paper [37] . 

The preceding analysis allows us to discuss the order parameter for the Mott 
transition. As it has been already discussed by different authors ( [3], [1], [33]), one 
of the main problems in characterizing the Mott transition is the absence of any 
obvious order parameter. For the Hubbard model, Castellani et al have identified it 
with the number of doubly occupied sites, by mapping the system onto the Blume- 
Emery-Griffits model, which describes the mixed phases of ^He He. Moreover, 
in the context of Dynamical Mean Field Theory, which maps the system onto a self- 
consistent impurity Anderson model, it has been argued that the order parameter 
is related to the effective hybridization of the medium, which captures the localiz- 
ing to delocalizing character of the transition. For the extended Hubbard model, 
other possibilities have been discussed in the literature. Let us consider here the 
(i-density wave order parameter, as defined by Wilczeck and Nayak in any number 
of dimensions: for a class of systems, it is assumed the existence of the non-trivial 
correlation function: 

(cU,.Cfc) = if{k) , (5.16) 
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where the vector Q has components of magnitude n/a in every direction and f{k) 
is a real function that changes sign under rotations of 7r/2 around any axis. This 
correlation function breaks a number of symmetries, namely time inversion T, 7c/2 
rotations, and lattice translations by one lattice unit, but preserves the product or 
the square of two of them. As it was pointed out in ref [3l] , this correlation function 
can be considered as arising from an order parameter which halves the size of the 
Brillouin zone and naturally induces the Mott transition. A quantum state realizing 
this correlation may be thought of as a d-wave spin singlet. In fact, the the following 
BCS-like wave function 

* = n^^'^T + cIt + ^feTCl+QT)Ki4i + VkiCllO) (5.17) 

k 

is energetically favorable in the extended Hubbard model. The relevance of this 
order parameter for the EFT is as follows: on the one hand, we expect parity to be 
broken in the EFT if the (i-density wave order parameter acquires a non- vanishing 
expectation value at (or near to) the Mott transition, as a consequence of the CPT 
theorem and charge conjugation symmetry. On the other hand, it is clear from its 
definition that the order parameter induces a short (or intermediate) scale charge 
ordering, so that, at the corresponding energy scales, a mass gap in the spectrum 
may develop, leading to an unnatural EFT. In this case, one is led to expect new 
physics to arise, leading to an EFT that would look very different from the initial 
one, implying that a search for a new theory should be undertaken, with new degrees 
of freedom to correctly describe the systeiifl. 

As it has been remarked in [15], long spin lifetimes in Si-Ga As samples appear 
at densities values near those that yield the Mott transition, which suggests that 
the latter might play a role in the anomalous spin behavior. At first glance, this 
appears far from obvious, because the Mott transition seems to be a phenomenon 
associated to the charge sector of the theory only. However, the EFT can give us 
some information given that we can compare the relative weight in the Hamiltonian 
of both the Dessehaus and the Coulomb terms at doping levels near criticality. 
Given that we know the scaling of both interaction terms, we can argue that both 
should be considered as equally important since at the strong fixed they become 
marginal simultaneously (since both scale to zero as when the energy scales down 
as s). This RG constraint connects both charge and spin sectors of the theory 
and, therefore, does not exclude the possibility that an spin effect, such as the one 
observed in the anomalous spin lifetimes, could be related to a Coulomb effect in 
the strong coupling regime. 

^ In fact, in two dimensions the d-density wave has been associated to the staggered flux phase 
[36] , which suggest the couphng of different Luttinger hquids, such that the main characteristics 
of each one (e.^.jthe gap formation and the Chiral symmetry breaking) are preserved. 
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6 Conclusions 



In this article wc have presented EFTs for electrons in crystalline systems and have 
studied the effect of quenched impurities (low disorder effects) in them. Based on 
experimental evidence and consistency and plausibility arguments, we have consid- 
ered the situation in which impurities form arrays in mesoscopic regions for which 
some or all of the discrete symmetries of the microscopic lattice are preserved in 
the long distance limit. We have also discussed EFTs for generic systems with sub- 
stitutional impurities, although we did not get into details as to how to deal with 
a general case. We beheve that the EFT point of view is a good complement to 
the one based on microscopic models for studying strongly correlated electron sys- 
tems in crystals, providing support to them and bringing in a useful bridge between 
strongly correlated problems and perturbative band theory calculations. We have 
obtained the EFT for 5'i-doped GaAs systems, obtaining from the symmetries of the 
effective fields all possible single particle terms in the effective Hamiltonian, which 
reproduce the known terms of the pseudo-potential calculations of band structure 
approaches and the so called Dyakonov-Perel Hamiltonian for the spin-orbit cou- 
pling. Remarkably, the resulting low-energy effective theory is described in terms of 
three (six chiral) one- dimensional Luttinger liquid systems and their corresponding 
intervalley transitions. 
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